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EVALUATION OF IMPORTANCE OF LATERAL ACCELERATION 
DERIVATIVES IN EXTRACTION OF LATERAL-DIRECTIONAL 
DERIVATIVES AT HIGH ANGLES OF ATTACK 

By Luat T. Nguyen 
Langley Research Center 

SUMMARY 


A theoretical investigation was conducted to determine the importance of the lateral 
acceleration (/S) derivatives in the extraction of lateral-directional stability derivatives 
for swept wing airplanes at high angles of attack. Representative values of lateral accel- 
eration derivatives in yaw and roll were used in a computer program to 

generate representative flight motions at several angles of attack and altitudes. The 
computer -generated motions were then subjected to a parameter identification process 
based on a modified Newton-Raphson method. Two identification techniques were evalu- 
ated, one which included the /3 derivatives and one which neglected them. The results 
of the identification procedures were compared with the actual values used to generate the 
flight motions. 


The results of the study indicate that omission of the /S derivatives from mathe- 
matical models used in derivative -extraction techniques can produce erroneous values for 
the lateral-directional stability derivatives particularly at high angles of attack, where the 
/3 derivatives are large. The largest errors occur in the dynamic derivatives, but large 
errors may also occur in the static derivatives for cases in which the ^ derivatives 
have large effects on the flight motions of the airplane. In addition, the resulting identi- 
fied mathematical models provide poor motion prediction as well as erroneous predictions 
of dynamic stability characteristics. These results strongly indicate that the effects of 
)3 derivatives shoifld be considered in any attempt to extract lateral -directional aerody- 
namic parameters at high angles of attack. 


INTRODUCTION 

The extraction of aerodynamic parameters from flight-test data has received 
increased attention over the past several years. Interest has been spurred primarily by 
the development of sophisticated parameter identification techniques which operate effi- 
ciently by using digital computers and which have been found to produce results superior 



to those obtained from earlier estimation methods. A detailed discussion of the advan- 
tages of the new techniques is given in reference 1. These new techniques are generally 
classified as output-error methods and include the so-called modified Newton-Raphson 
and quasi -linearization methods. (See refs. 2 and 3.) Recently, a great deal of interest 
has been generated toward applying these parameter identification techniques to obtain 
estimates of aerodynamic derivatives for fighter airplanes in the high-angle -of -attack 
flight regime. 

A primary consideration in this application is the mathematical model used in the 
identification process. It is well known that output-error methods can give very ques- 
tionable results in the presence of process noise such as modeling errors. Up to the 
present time, the lateral -directional stability derivatives due to lateral acceleration (/s) 
have not been included in mathematical models used in the estimation of lateral -directional 
aerodynamic parameters. Such derivatives have usually been omitted because they are 
generally small at low angles of attack; moreover they are not easily estimated or meas- 
ured in the wind tunnel. As indicated in reference 4, past wind-tunnel investigations 
employing testing equipment designed specifically to obtain the lateral acceleration deriv- 
atives have shown that large values for the derivatives C, . and C„ • can be obtained 

for configurations with highly swept wings at high angles of attack. The past studies have 
also shown that these derivatives can have a major effect on the dynamic lateral- 
directional characteristics of a typical fighter at high angles of attack. 


The present investigation was conducted to examine the importance of the lateral 
acceleration derivatives in the lateral -directional stability derivative identification 
process at high angles of attack. The effect of omitting these derivatives on the resulting 
parameter estimates was investigated by use of two linear dynamic models. One of the 


models included the effects of C, . and 

h 


whereas the second model represented the 


conventional approach in which these derivatives are neglected. The estimated aerody- 


namic parameters which resulted with and without the C, . and Cn • terms were com- 
pared with the actual values of the respective parameters. 


SYMBOLS 


The lateral-directional derivatives presented herein are referred to the body axis 
system shown in figure 1. Dimensional quantities are presented in both U.S. Customary 
Units and the International System of Units. The subscript o indicates nominal value, 
and a dot over a symbol denotes a time derivative. The primed dimensional derivatives 
^for example, Lp, N^, Njg, etc.^ are defined by equations (A9) to (A24) in the appendix. 
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wing span, m. (ft) 
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Cn 
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Cy 
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lift coefficient, Aerodynamic j ift for ce 

qS 


rolling-moment coefficient about X body axis. 


Aerodynamic rolling moment 


pitching-moment coefficient about Y body axis. 
Aerodynamic pitching moment 
qSb 

yawing-moment coefficient about Z body axis, 


Aerodynamic yawing moment 


X-axis force coefficient along positive X body axis. 
Aerodynamic X-axis force 
qS 


Y-axis force coefficient along positive Y body axis. 
Aerodynamic Y-axis force 

qS 


Z-axis force coefficient along positive Z body axis. 

Aerodynamic Z-axis force 
qS 

wing mean aerodynamic chord, m (ft) 
acceleration due to gravity, m/sec^ (ft/sec^) 
altitude, m (ft) 

moment of inertia about X body axis, kg-m^ (slug-ft^) 

product of inertia with respect to X and Z body axes, kg-m^ (slug-ft^) 

moment of inertia about Y body axis, kg-m^ (slug-ft^) 
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moment of inertia about Z body axis, kg-m^ (slug-ft^) 
airplane mass, kg (slugs) 
period, sec 

airplane roll rate about X body axis, ra.d / sec or deg/ sec 

airplane pitch rate about Y body axis, rad/sec or deg/sec 

free-stream dynamic pressure, N/m^ (Ib/ft^) 

yaw rate about Z body axis, rad/ sec or deg/sec 

wing area, m^ (ft^) 

time, sec 

time to damp to one -half, sec 

component of airplane velocity along X body axis, m/sec (ft/sec) 
airplane resultant velocity, m/sec (ft/sec) 

component of airplane velocity along Y body axis, m/sec (ft/sec) 
component of airplane velocity along Z body axis, m/sec (ft/sec) 
airplane orthogonal body axes (see fig. 1) 
angle of attack, rad or deg 
angle of sideslip, rad or deg 

aileron deflection, positive for left roll, rad or deg 
elevator deflection, positive for nose -down control, rad or deg 


rudder deflection, positive for left yaw, rad or deg 





Euler angles, rad or deg 
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METHOD OF ANALYSIS 


A schematic of the overall approach used in the investigation is shown in figure 2. 
Initially, mass and aerodynamic characteristics representative of a swept -wing fighter 
configuration were used as inputs to a nonlinear, stx-degree-of -freedom computer pro- 
gram which generated time histories representative of flight data. The simulated flight 
data were then used in a Newton-Raphson parameter identification program which identi- 
fied unknown lateral -directional stability derivatives by use of linear state-space models 
of the airplane. Two mathematical models were used for the parameter identification 
process. The first model properly represented and whereas the second 

model omitted these terms. The output of the Newton-Raphson program was a set of 
estimated aerodynamic derivatives extracted from the simulated flight data. Since the 
exact values for the aerodynamic derivatives were known inputs, the effects of mathemat- 
ical model variations could be easily examined. 
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Simulation of Flight Data 

The airplane configuration used in the investigation was the represehtative swept- 
wing fighter used in reference 4. The mass and dimensional characteristics of the con- 
figuration are listed in table I, and the aerodynamic data are presented in table II. The 
aerodynamic data were based on static and dynamic wind-tunnel tests of several configu- 
rations.. It should be noted that the derivatives given in reference 4 were converted from 
the stability axis system to the body axis system for the present study. Furthermore, the 
study reported in reference 4 was restricted to lateral-directional considerations. In 
order to represent flight motions of an airplane with six degrees of freedom, represen- 
tative values of longitudinal stability derivatives were chosen for the present study. 

The initial values of Cj^ were identical to those used in reference 4. Representative 
values of the control derivatives were also used in the investigation. 


The variations of the lateral-directional stability derivatives with angle of attack 

are shown in figure 3. The derivatives Cj^. and . were relatively small at 

^ /3 

a = 10°, but the magnitudes of the derivatives increased rapidly as a was increased. 


At higher angles of attack, two variations of 


Cn and 
°/3 


'Ir. 


were used to represent a 


range of static stability characteristics noted for fighter -type airplanes. The first vari- 
ation (labeled case (a)) was representative of cortfigurations which lose directional sta- 
bility and effective dihedral at high angles of attack. This type of configuration usually 
exhibits a directional divergence ("nose slice") near stall. (See ref. 5.) The second var- 
iation (labeled case (b)) was representative of configurations which maintain stability at 
high angles of attack and usually have good stall characteristics. 

The flight conditions for which data were generated are presented in table IE. The 
initial conditions involved steady, level flight at angles of attack of 10°, 20°, and 30° at 
altitudes of sea level, 7620 m (25 000 ft), and 15 240 m (50 000 ft). 

The dynamic lateral -directional stability characteristics of the airplane at these 
flight conditions are given in table IV. Part (a) of table IV presents values obtained when 
the terms and were included in the equations of motion to obtain the period 

and damping of the various modes. For comparison, stability characteristics obtained 
when both and were assumed to be zero are listed in part (b) of table IV. The 

differences between the calculated stability characteristics for the two sets of data are 
small at a = 10° ^where the values of and are smallj , fairly large at 

a = 20° (where the derivatives are larger), and very large at o: = 30° (where the 
/3 derivatives were very large). These results show that the lateral acceleration deriv- 
atives can have a dominant effect on dynamic lateral-directional stability, particularly at 
high angles of attack, where these derivatives are large. 


7 



It will be noted that the results of table IV for a = 30° (for example, runs 3(a) 
and 3(b)) show markedly different stability characteristics for the two sets of 

and C, used. For those designated (a), the unstable stability derivatives result in an 

^i3 

aperiodic divergence ("nose slice"). The results obtained for those designated (b) indi- 
cate a stable airplane at a = 30°. 

The aerodynamic and mass characteristics were used as inputs for a six-degree - 
of -freedom (6 DOF) nonlinear computer program which generated time histories of 
motions following control inputs. At each flight condition lateral -directional motion was 
initiated with the application of a rudder doublet followed by an aileron doublet. The air- 
plane was then allowed a period of free motion during which the lateral -directional con- 
trols were neutralized. In generating these motions, an attempt was made to excite the 
lateral-directional modes sufficiently for proper derivative identification, whereas the 
motions were restricted so that the small perturbation assumption of the mathematical 
model was not invalidated. The motion program integrated the equations of motion at a 
fixed interval of 0.01 sec, and selected lateral-directional flight variables were provided 
as output on punched cards at 0.04 sec intervals for input to the identification program. 
Time histories of a typical flight are given in figure 4. The data show very little longi- 
tudinal motion; at the same time the lateral-directional modes were sufficiently excited 
without the motion becoming large enough to invalidate the linear, small perturbation 
model. 


Estimation of Aerodynamic Parameters 

The technique used for estimating the aerodynamic parameters from the computer- 
generated flight data was a modified Newton-Raphson method. A detailed description of 
the technique and computer program is found in reference 2. The technique computes a 
maximum -likelihood estimate of the unknown parameters of a linear, state-space model. 
Because measurement noise was not considered in this investigation, the maximum- 
likelihood estimate obtained is equivalent to one which minimizes the mean-square differ- 
ence between the computed and measured response. The flight parameters that were 
input into the estimation program were the response variables p,r,/3,0 and the control 
deflections 6g^ and 6j.. These parameters are commonly measured during flight 
testing and are normally used for identification of aircraft lateral -directional aerody- 
namic parameters. No measurement of ^ was assumed to be available. Because the 
measurements were ideal (that is, no measurement noise was included in the generation 
of the flight data), differences between the measured and computed response variables 
p,r,/3,0 were all weighted equally in the minimization process. 
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The linearized equations of motion used in the identification process are described 
in detail in the appendix. In the development of the mathematical model, the classical 
lateral -directional equations, including terms which were the result of the ^ derivatives, 
were put into state-space form to be compatible with the identification program. In this 
form, all the aerodynamic parameters appear as coefficients of the state variables 
p,r,/3,0 and of the controls 6^,6^,. In matrix form, the equations are as follows: 




^/3 

% 

n; 



y; 


1 

tan 6 ^ 

0 


0 


P 


La 

^6 

L>r 

n’w, 

0 


r 

+ 


NSr 

0 


|3 



Y'sr 

0 


0 


0 

0 


With this model, the |3 derivatives appear in all the primed derivatives but are explicitly 

identifiable as components of the coefficients of the bank angle 0. For example, the 

lV and N', terms are related to these derivatives as 

0 0 


and 




1 


j ^XZ 

_\Vo 

Vz 

1 

_/_g_ 

f2 

1 ^XZ 

Vo 


cos 6 




cos 9 . 






With the |3 derivatives identified, the other derivatives can be determined from the 
remaining primed derivatives by using the relationships given in the appendix. 

Two models were used in the identification process. The first model included the 
values of C^ . and estimates were obtained for the following parameters; 

Lp’ l;, L'^, l’^, n;^, n;, N’^, n’^, y; 3 , L'g^, N’^,and N’^. The values 

of the derivatives Cv , Cv , and were assumed to be 0, and C-rr was known. 

^ Yp’ Y5^ 

Thus, Yp, yJ^, Yg , Yg , and Y*, were also known and hence not included in the set 

3, I* r 

of parameters to be identified. Once estimates of the dimensional primed derivatives 
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were obtained, the conventional nondimen sional stability derivatives were derived by 
using the relationships given in the appendix and symbols section. 

The second model used for parameter identification was one in which the values of 
C, . and C„. were assumed to be 0. The values of L*. and N*. were therefore 

assumed to be 0, and the computer-generated flight motions (with effects of the /3 deriv- 
atives present) were subjected to the conventional derivative identification process in 
which the /3 terms are neglected. 

All the estimation runs were 10 sec in duration with a sample rate of 0.04 sec, for 
a total of 250 data points. In each case, the estimation program was allowed to operate 
until a convergence to a solution was established. Convergence always occurred within 
the first 10 iterations; however, to assure complete convergence, at least 15 iterations 
were made in every case. At the point at which the program was stopped, values of the 
gradient of the minimization fimction were typically on the order of 10"10 or less and 
indicated an accurate convergence. 

RESULTS AND DISCUSSION 

Effect of /3 Derivatives on Standard Deviation and Data Fit 

Presented in table V are values of the primed dimensional lateral -directional 
derivatives and calculated estimates of the standard deviation for the two extraction 
techniques. The standard deviation is a measure of uncertainty in the extracted param- 
eter, and a small value of the deviation in comparison with the magnitude of the parameter 
indicates a high degree of certainty in the estimate. The data of table V indicate that the 
standard deviations for the extracted derivatives were extremely small for all flight con- 
ditions for both identification techniques. Also, figures 5 to 8 show that the motions gen- 
erated when either set of derivatives is used agree very closely with the actual "flight" 
motions from which the derivatives were extracted. 

For the case of the unmodeled ^ derivatives, the foregoing results are completely 
misleading. A high degree of certainty is indicated by the standard deviation and time- 
history fits, even though omission of the ^ derivatives caused considerable errors in 
their estimated values. Moreover, when these derivatives were used to compute dynamic 
characteristics or to predict motions other than those used to extract the derivatives 
themselves, the results obtained were in error. (See discussion in the following section.) 

Comparison of Extracted Derivatives With Actual Values 

In figures 9 to 11, the values of the conventional nondimens ional lateral-directional 
stability derivatives obtained by the two identification techniques are compared with the 
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actual values used in the computer -generated motions . The comparisons are shown at 
a = 10°, 20°, and 30° for altitudes of sea level, 7620 m (25 000 ft), and 15 240 m 
(50 000 ft). In general, the comparisons show that the errors caused by neglecting the 
)3 derivatives were small at o! = 10° (where the /3 derivatives were relatively small), 
larger at a = 20°, and very large at a = 30° (where the /3 derivatives were large). 


The largest errors were obtained in the dynamic derivatives, as might be expected. 
In addition, large errors in the static derivatives were also obtained for cases where the 
/3 derivatives were large. For example, for a = 30° at sea level (fig. 9), omission of 
the ^ derivatives resulted in extracted values of and Cj^ which were markedly 

different from the actual values. In fact, for case (a) ^unstable values of and 

at Of = 30°^ the estimates were opposite in sign to the actual values. This particular 
result was caused possibly by the fact that the /3 derivatives provided substantial 
damping of the lateral-directional modes at a = 30°. (See table IV.) 


When the flight motion created with the jS derivatives was subjected to the identi- 
fication process which omitted these terms, the technique forced the values of C, 

^/3 

and Cnp to become stable in order to provide the increase in stability previously pro- 
duced by the /3 derivatives. The values and trends estimated by the technique which 
neglected • and . were markedly different from those for the actual data used in 

generating the flight motions and indicated thatt severe errors are produced by omitting 
these terms in extraction processes at high angles of attack. 


Also presented in figure 9 are values indicating the combined dynamic lateral- 
directional parameters normally measured during wind-tunnel forced-oscillation tests 
in roll and yaw. (See ref. 6.) The various combined derivatives ^such as 
C^r " obtained rather than pure derivatives because of the physical con- 

straints of wind-tunnel test procedures and equipment. It is interesting to note that the 
derivatives extracted by the identification technique which omitted the ^ derivatives 
agree reasonably well with these values. The reason for this agreement can be seen by 
examination of the equations relating the primed derivatives to the dimensional derivatives 
as shown in the appendix (eqs. (A13) to (A16)). For example, is given by 



^XZ 

ix 


N^^-cos a 


o + YJ 
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Because Ixz/^X normally small, the dominant term in the right-hand side 

of the equation is L^. - cos a^. Inspection of the remaining equations led to similar 
conclusions; that is, 


Lp + sin 


n: 


p 


Np + sin Q-o 


Nr « Nj, - cos o-^ 

These results indicate that a derivative -extraction process which omits the 
/3 derivatives at high angles of attack (where they are large) produces values of the 
dynamic derivatives which tend to agree, at least in trend, with the combined derivatives 
obtained in wind-tunnel forced-oscillation tests. As stated earlier, however, such an 
extraction process not only provides erroneous estimates of the stability derivatives but 
also poor predictions of the airplane dynamic stability characteristics as defined by the 
period and time to damp to one-half amplitude of the characteristic modes. Shown in 
table VI are comparisons of lateral-directional stability characteristics computed by 
using the actual aerodynamic data and computed by using extracted data from the non- 
/3-modeling technique. Note that the predicted modal characteristics differ considerably 
from the actual characteristics for the cases in which a = 30°, where the |3 derivatives 
are very large. 

The fact that neither the modal characteristics nor the aerodynamic derivatives 
were correctly obtained when the /3 derivatives were omitted from the extraction pro- 
cedure indicates that although the resulting identified mathematical models generate 
motions which closely match the flight motions from which they were extracted, these 
models could not be used to predict motions for other flight conditions or for different 
control inputs. For example, figure 12 presents comparative time histories of responses 
to an aileron doublet applied at the conditions of test run 3(a) {a = 30°, h = 0). A set of 
responses computed by using the results obtained from the non-/3-modeling estimation 
technique is compared with the actual responses. As can be seen, the predicted responses 
do not agree well with the actual responses. 

These resvilts indicate that the omission of /3 derivatives from the derivative- 
extraction procedure at high angles of attack, where these derivatives are large, results 
in identified mathematical models which provide poor motion prediction as well as 
erroneous predictions of dynamic mode characteristics. 
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Also plotted in figures 9 to 11 are the results obtained when the derivatives 
were included in the identification model in terms of L*^ and N'^. As expected, the 
results agree closely with the actual values. It sho^^ld be noted that the primary effects 
of the /3 derivatives appear in terms of L*^ rather than N^. This fact is reflected 
in the small values of n', as shown in table V: note that in some cases nV was fixed 
at 0 and not identified. On the other hand, values of L'^ were large enough to be always 
easily identifiable. As shown in the figures, the values of 0^^ and C^. calculated 

from the l', and N*. estimates were in excellent agreement with the actual values. 




<P 


The data also reveal that the estimates of the control derivatives agree reasonably 
well with the actual values for both identification techniques. This result might be 
expected since the nonmodeling of /3 derivatives should not have a major effect on the 
identification of control parameters. 


The data of figures 10 and 11 show that the extraction errors previously discussed 
for flight at sea level were present at higher altitudes, but the errors for and C, 

were less at the higher altitudes. The smaller error was probably caused by the fact 
that the /3 derivatives had a reduced effect on the stability of the airplane at higher 
altitudes. (See table IV.) 


CONCLUSIONS 

The results of a theoretical study to determine the importance of the lateral accel- 
eration derivatives in the extraction of lateral -directional stability derivatives at 
high angles of attack has produced the following conclusions: 

1. Omission of the /3 derivatives from mathematical models used to represent 
airplane flight conditions for derivative extraction at high angles of attack produces 
erroneous values for the lateral -directional stability derivatives at conditions where the 
/3 derivatives are large. 

2. The largest errors occur in the dynamic derivatives, but large errors may also 
occur in the static derivatives in cases where the /3 derivatives have large effects on 
the flight motions of the airplane. 

3. At high angles of attack, where the jS effects are large, the mathematical 
models identified when the /3 derivatives are omitted provide poor motion prediction as 
well as erroneous predictions of dynamic mode characteristics. 
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4. The results strongly indicate that the effect of /3 derivatives should be con- 
sidered in any attempt to extract lateral -directional aerodynamic parameters at high 
angles of attack. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 17, 1974. 
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APPENDIX 

LINEARIZED EQUATIONS OF MOTION USED 
IN IDENTIFICATION PROCESS 

The general three-degree-of -freedom linearized lateral -directional equations of 
motion are (see ref. 7, for example): 


/3 + r cos Q!q - P sin cos e^cp = Y^/3 + + YpP + Y^r + Yg 5^ + Yg^6^ 


(Al) 


. ^Z *■ ^Y 
p + — — qo>^ 


xz 

X 


(v 



= L^/3 + + LpP + L^r + Lg^6^ + Lg^6j, 


(A2) 


r + 



%P + 


^Z 





N^/3 + NpP + N^r + N 


0^ + 


N« 6, 


(A3) 


0 = p + tan 0^r + tan 6^(p 


(A4) 


The /3, p, and r equations can also be expressed as 


•6^ 


g 

Y. sin + Y„ -cos + Y„ Yg 

/3= ^,3 + 2 Pp + 2 £r + -5 0 +_£a^5a. r 

1-Y^ 1-Y^ 1-Y^ 1-Y^ 1-Y^ 1-Y^ 


6. 


(A5) 
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APPENDIX - Continued 





16 



APPENDIX - Continued 






(A7) 
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APPENDIX - Continued 


By using matrix notation, equations (A4), (A5), (A6), and (A7) can be written as 



l: 


n: 


y' 

r 

tan 0Q 



l’ 


P 



LqI 

Or 

< 


r 

+ 

N'Sa 

N6r 






YSr 

% tan 0Q 




0 

0 


(A8) 


For this investigation = 0, and the flight conditions were such that Qq = 0. As a 
result, the primed derivatives can be written in terms of the dimensional derivatives as 
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APPENDIX - Continued 


N* = ^ 

t 2 


"XZ 




n 


^xz 


+ Np] + ^ L^(sin + Yp) + N^(sin o-^ + Yp) 


(A14) 


t' _ 1 

:5 


1 - 


^XZ 




(^r ^ '^rj + L^(-cos <»„ Yr) + ^ N^(-cos o,^ + Y,) 


(A15) 


n’ = L_ 


^XZ 


h^^z 


+ ”rj + ^ (-“" “o ^r) (-“= “o + ^r) 


(A16) 


La = — V 

t 2 


1 - 


^XZ 




|Lg H-^Ng | + (LAYg +^^NaY5 
^ Ix °al \ p 5a Ix ^ aj 


(A17) 


Nk = 


a j2 
^XZ 


-.Ng W(^^L-Y, .N.Y, 

Iz 5al 1 Iz ^ ^a ^ ^ai 


^X^Z 


(A18) 


La = - ^ 


5 r ^ 


1 - 


XZ 

^Z 




(A19) 


Nk = 


^6 ~ 9 

r l4 


1 - 


XZ 

%^Z 




(A20) 


Yp = sin Q!^ + Yp 


(A21) 
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APPENDIX - Concluded 


Yj. = -cos o-Q + Yj. 


(A22) 




(A23) 


y , = cos e_ 
Vo 


(A24) 


Note that the lateral acceleration derivatives appear in all the primed roll and yaw 
derivatives and that the parameters and n’^ are composed of them exclusively 

Thus, identifying and N*^ makes possible the determination of . and C 

directly. ^ 
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TABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS OF AIRPLANE 


Weight, N (lb) 

Moments of inertia, kg-m^ (slug-ft^) 



^XZ 


Wing dimensions: 

Span, m (ft) 

Area, m^ (ft^) 

Mean aerodynamic chord, m (ft) 


101 641.4 (22 850 ) 

18 915.3 (13 951 . 39 ) 
94 906.0 (70 000 . 0 ) 
117 940.4 (86 989 . 51 ) 
2071.0 ( 1527 . 50 ) 


11.61 ( 38 . 1 ) 
61.50 ( 662 . 0 ) 
5.28 ( 17 . 32 ) 
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TABLE II. - AERODYNAMIC COEFFICIENTS 


Designations (a) and (b) on 



and 



indicate the two sets of data used for runs 3, 6, and 9 


Coefficient 

Values of a of - 

0° 

10° 

20° 

30° 

40° 

Longitudinal 


-0.05 

-0.03887 

-0.13045 

-0.31407 

-0.46031 


0 

-.41303 

-.89882 

-1.3360 

-1.6917 

Cm 

.04 

0 

-.04 

-.08 

-.12 

Cx 

\ 

0 

0 

0 

0 

0 

C y 

0 

0 

0 

0 

0 


-15.424 

-17.967 

-21.078 

-15.155 

-7.707 


0 

0 

0 

0 

0 

Cz 

0 

0 

0 

0 

0 

■"»e 

-.00977 

-.01494 

-.01505 

-.01637 

-.01384 

Lateral directional 

C, (a) 
^/3 

-0.0008433 

-0.001159 

-0.001474 

0.00050 

0.002474 

C, (b) 
^/3 

-.0008433 

-.001159 

-.001474 

-.001366 

-.001258 


.001095 

.0008112 

.0005278 

-.000866 

-.00226 

c„/f) 

.001095 

.0008112 

.00,05278 

.000366 

.0002042 

Cy 

/3 

-.009948 

-.009948 

-.009948 

-.004992 

-.000035 

Cy • 
^13 

0 

0 

0 

0 

0 
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TABLE II.- AERODYNAMIC COEFFICIENTS - Concluded 


Coefficient 

Values of O' of - 

0° 

10° 

20° 

30° 

40° 

Lateral directional - Concluded 

c,. 

0.18876 

-0.15467 

-0.49810 

-1.1062 

-1.7143 

Cn- 

-.02614 

.01334 

.05282 

.51603 

.97924 

Cy 

0 

0 

0 

0 

0 


-.23229 

-.17801 

-.12373 

. 

.08990 

.30353 

Cn 

"P 

.02147 

.002115 

-.01724 

-.09036 

-.16348 


0 

0 

0 

0 

0 

6r 

.17648 

.10212 

.027761 

.009641 

-.00848 

Cn 

Up 

-.13773 

-.1720 

-.20627 

-.20990 

-.21353 

Cy. 

0 

0 

0 

0 

0 

C/ 

-.00138 

-.00102 

-.00067 

-.00045 

-.0002 

Cn. 

^a 

0 

0 

0 

0 

0 

Cv 

Y6p 

.00369 

.00399 

.00344 

.00173 

.00075 

C7 

^6r 

0 

0 

0 

0 

0 

S 

Op 

-.00184 

-.00202 

i 

-.00167 

-.00069 

-.00021 
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st runs 3, 6, and 9 correspond 

and C„ shown in table II 

iS ^/3 


(ft) 

V, m/sec (ft/sec) 

(0) 

80.19 (263.09) 

(0) 

53.10 (174.21) 

(0) 

41.81 (137.18) 

(0) 

41.81 (137.18) 

(25 000) 

119.73 (392.80) 

(25 000) 

79.28 (260.10) 

(25 000) 

62.43 (204.81) 

(25 000) 

62.43 (204.81) 

(50 000) 

209.94 (672.37) 

(50 000) 

135.71 (445.23) 

(50 000) 

106.86 (350.58) 

(50 000) 

106.86 (350.58) 
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TABLE IV.- AIRPLANE DYNAMIC LATERAL-DIRECTIONAL STABILITY 

CHARACTERISTICS AT STUDY FLIGHT CONDITIONS 

^Designations (a) and (b) on test runs 3, 6, and 9 correspond 
to the two data sets for C^^ and C^^ shown in table IlJ 


(a) 4 0; C^^ 4 0 


Test 

Aperiodic modes 

Oscillatory mode 

rxin 

tj^/ 2 , sec 

H/2’ sec 

P, sec 

ti/ 2 , sec 

1 

14.91 

0.429 

4.538 

1.409 

2 

5.393 

.915 

5.403 

.8907 

3(a) 

-3.821 

.254 

53.93 

2.369 

3(b) 

6.754 

1.503 

9.634 

.5401 

4 

21.91 

.72 

4.21 

1.79 

5 

7.97 

1.44 

4.74 

1.30 

6(a) 

-2.58 

.34 

49.13 

4.52 

6(b) 

9.83 

1.98 

6.06 

.83 

7 

37.17 

1.317 

3.995 

2.821 

8 

13.56 

2.515 

4.476 

2.204 

9(a) 

-1.55 

.451 

78.30 

10.36 

9(b) 

16.62 

3.243 

5.218 

1.430 


(b) 

Cr, . = 0; C, . = 

0 


Test 

Aperiodic modes 

Oscillatory mode 

run 

: 1 
H/2’ 

ti/ 2 , sec 

P , sec 

H/2’ 

1 

14.73 

0.451 

4.345 

1.777 

2 

5.282 

.956 

4.479 

3.714 

3(a) 

-.569 

.979 

30.65 

5.69 

3(b) 

6.50 

1.49 

4.69 

-1.66 

4 

21.78 

.73 

4.12 

2.29 

5 

7.89 

1.46 

4.41 

5.34 

6 (a) 

-.63 

.92 

46.36 

7.75 

6(b) 

9.67 

2.05 

4.78 

-2.36 

7 

37.1 

1.33 

3.97 

3.61 

8 

13.5 

2.52 

4.38 

8.96 

9(a) 

.857 

-.678 

80.4 

12.7 

9(b) 

16.5 

3.31 

4.85 

-3.91 
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TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS 


(a) a = 10*^; h = 0 


Parameter 

Unmodeled 

jS effects 

Modeled (. 

effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-2.3673 

0.0009 

-2.2726 

0.0007 


2.7979 

.0004 

2.8093 

.0004 

^/3 

-10.0707 

.0001 

-9.5129 

.0001 




-.1792 

.0006 

L’e 

-9.8037 

.0003 

-9.6780 

.0003 

^6 

-.3675 

.0007 

-.3447 

.0007 

Np 

-.0752 

.0018 

-.0285 

.0016 

Nr 

-.3270 

.0018 

-.2634 

.0010 


.6789 

.0004 

.9530 

.0003 




.0018 

.0046 

Ng 

°a 

-.2459 

.0007 

-.1547 

.0005 

Ng 

Or 

-2.7796 

.0012 

-2.7645 

.0011 

y' 

P 1 

-.2108 

.0015 

-.1689 

.0009 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 

(b) a = 20*^; h = 0 


Parameter 

Unmodeled 

/3 effects 

Modeled ^ effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-2.2062 

0.0008 

-2.1162 

0.0008 

4 

3.1723 

.0003 

3.5756 

.0003 


-6.4493 

.0002 

-5.0728 

■ 

.0001 




-.6306 

.0007 


-2.5916 

.0010 

-2.5133 

.0010 

4r 

-.2994 

.0012 

-.2245 

.0013 


-.0748 

.0018 

1 

-.0359 

.0018 

4 

-.3617 

.0009 1 

-.2286 

.0007 

% 

.0569 

.0007 

.2157 

.0008 

N'w, 



.0019 

.0045 

Ng 

®a 

-.0273 

.0021 

-.0441 

.0020 

% 

Uj. 

-.9885 

.0019 

-1.0002 

.0017 

y' 

^)3 

-.3062 

.0009 

-.1076 

.0007 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(c) a = 30°; case (a), h = 0 


Parameter 

Unmodeled 

^ effects 

Modeled (. 

effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 


-2.5012 

0.0004 

-2.6800 

0.0004 

4 

4..2980 

.0002 

5.5711 

.0001 


-2.7270 

.0001 

1.4812 

.0001 




-1.3440 

.0003 


-.9977 

.0015 

-1.0883 

.0013 

LOr 

-.2539 

.0023 

-.1059 

.0017 

Np 

.0671 

.0012 

.0948 

.0012 

Nr 

-.5092 

.0005 

-.4753 

.0004 

n’^ 

.1166 

.0006 

-.2930 

.0004 

n’. 



.0680 

.0012 


-.0561 

.0028 

-.0222 

.0027 

N’Cr 

-.2560 

.0045 

-.2514 

.0038 

Y’ 

P 

-.0707 

.0006 

-.0301 

.0007 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(d) a = 30°; case (b); h = 0 


Parameter 

Unmodeled 

/3 effects 

Modeled ^ effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

■ 

Lp 

-2.3405 

0.0004 

-2.6388 

0.0004 

4 

4.1234 

.0002 

5.5176 

.0001 


-5.6868 

.0002 

-2.9549 

.0001 




-1.2843 

.0005 


-1.0005 

.0015 

-1.0628 

.0013 

4r 

-.1750 

.0031 

-.1146 

.0027 

Np 

-.0585 

.0014 

.1086 

.0009 

4 

-.3651 

.0006 

-.5053 

.0004 


.2263 

.0005 

.0678 

.0007 


^ — 


.0751 

.0018 


-.0790 

.0026 

-.0175 

.0022 

NSr 

-.2346 

.0050 

-.2518 

.0045 

Y’ 

-.3888 

.0005 

-.0431 

.0004 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(e) a = 10^^; h = 7620 m (25 000 ft) 


Parameter 

Unmodeled 

^ effects 

Modeled ^ effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 


-1.4971 

0.0011 

-1.4993 

0.0011 

l; 

1.9585 

.0003 

1.8433 

.0003 


-9.6108 

.0001 

-9.6758 

.0001 




-.0902 

.0008 

^6a 

-9.4561 

.0004 

-9.3159 

.0003 


-.3108 

.0004 

-.3098 

.0004 

Np 

-.0168 

.0030 

-.0219 

.0029 

Nt 

-.1801 

.0007 

-.1788 

.0007 


.9676 

.0004 

.9364 

.0004 

0 



0 (fixed) 


Ng 

-.1416 

.0009 

-.1548 

.0009 

^6 

Or 

-2.7576 

.0007 

-2.7569 

.0007 

Y’ 

^/3 

-.1108 

.0009 

-.1133 

.0009 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 

(f) a = 20°; h = 7620 m (25 000 ft) 


Parameter 

Unmodeled 

^ effects 

Modeled ^ effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 


-1.4468 

0.0014 

-1.4148 

0.0013 

Lr 

2.2520 

.0003 

2.3769 

.0003 


-5.8514 

.0003 

-5.3193 

.0002 




-.2804 

.0006 


-2.5535 

.0014 

-2.5149 

.0014 


-.2092 

.0010 

-.1569 

.0010 


-.0373 

.0027 

.0001 

.0180 

Nr 

-.1979 

.0009 

-.0238 

.0026 


.2029 

.0010 

-.1542 

.0008 




.2184 

.0010 

Ng 

°a 

-.0306 

.0026 

-.0441 

.0025 

N’s 

Or 

-.9822 

.0015 

-1.0015 

.0013 

y’ 

-.1080 

.0013 

-.0731 

.0009 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(g) a = 30*^; case (a); h = 7620 m (25 000 ft) 



Unmodeled 

^ effects 

Modeled ^ effects 

Parameter 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-1.8445 

0.0005 

-1.7521 

0.0005 

Lr 

3.0631 

.0003 

3.6494 

.0002 


1.5987 

.0003 

1.2553 

.0003 

t/ 



-.5728 

.0011 

<t> 



L's 

®a 

-1.4961 

.0010 

-1.0478 

.0010 

L6r 

-.1425 

.0030 

-.0653 

.0027 

% 

.2314 

.0011 

.0713 

.0014 

Nr 

-.7828 

.0006 

-.3350 

.0005 


-.6284 

.0007 

-.3079 

.0347 

.0011 

.0025 



Ns 

.1185 

.0017 

-.0197 

.0022 

N’sr 

-.3047 

.0056 

-.2544 

.0048 

Y* 

^/3 

.0459 

.0015 

-.0303 

.0016 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(h) a = 30°; case (b); h = 7620 m (25 000 ft) 


' 

Parameter 

Unmodeled 

effects 

Modeled 

/3 effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-1.6048 

0.0006 

-1.7566 

0.0006 

Lr 

3.1503 

.0003 

3.6604 

.0002 


-3.8954 

.0003 

-3.0854 

.0002 




-.5719 

.0009 

L’s 

°a 

-1.0242 

.0018 

-1.0545 

.0017 


-.1043 

.0040 

-.0653 

.0038 

Np 

-.0540 

.0022 

.0723 

.0013 

Nr 

-.2150 

.0009 

-.3355 

.0006 

Njs 

.1365 

.0008 

.0754 

.0008 

< 



.0337 

.0042 


-.0488 

.0028 

-.0177 

.0023 

n’s 

Or 

-.2327 

.0060 

-.2544 

.0051 

y' 

^/3 

-.3266 

.0006 

-.0278 

.0005 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(i) a = 10°; h = 15 240 m (50 000 ft) 


Parameter 

^ 

Unmodeled 

^ effects 

■ 

Modeled 

/3 effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-0.8720 

0.0018 

-0.8708 

0.0018 

Lr 

1.0157 

.0004 

1.1068 

.0004 


-9.9465 

.0002 

-9.6943 

.0002 




-.0313 

.0025 

L 5 

°a 

-9.5132 

.0004 

-9.5374 

.0004 


-.3477 

.0005 

-.2782 

.0005 

Np 

-.0107 

.0083 

-.0134 

.0080 

4 

-.1351 

.0016 

-.0981 

.0012 

N/3 

.9089 

.0009 

.9330 

.0008 

n’ 



0 (fixed) 


N 5 

®a 

-.1213 

.0013 

-.1318 

.0013 

4 

r 

-2.7660 

.0009 

-2.7455 

.0009 

Y’ 

^/3 

-.0695 

.0011 

-.0717 

.0011 
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TABLE V,- ESTIMATES OP AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 

(j) a = 20°; h = 15 240 m (50 000 ft) 


Parameter 

Unmodeled 

^ effects 

Modeled /3 effects 

Estimated 

value 

Deviation 

Estimated . 
value 

Deviation 

Lp 

-0.8218 

0.0022 

-0.8219 

0.0022 

Lr 

1.2189 

.0006 

1.4085 

.0005 

^13 

-5.7610 

.0004 

-5.3652 

.0003 




-.0992 

.0018 


-2.4817 

.0015 

-2.4947 

.0015 

Lg 

Or 

-.1540 

.0017 

-.1210 

.0017 

Np 

-.0237 

.0045 

-.0163 

.0048 

Nr 

-.1132 

.0021 

-.0766 

.0021 


.1538 

.0012 

.2206 

.0010 

•'V 



0 (fixed) 



-.0279 

.0038 

-.0423 

.0036 

NSr 

-.9843 

.0025 

-.9938 

.0026 

y’ 

-.1036 

.0013 

i 

-.0518 

1 

.0013 
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TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Continued 


(k) a = 30°; case (a); h = 15 240 m (50 000 ft) 



Unmodeled 

^ effects 

Modeled 

/3 effects 

Parameter 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 


-0.9612 

0.0012 

-1.0491 

0.0012 

Lr 

1.8218 

.0006 

2.1342 

.0005 


.6438 

.0003 

1.1889 

.0004 

0 

L'da 

-1.1254 

.0014 

-.2059 

-1.0803 

.0014 

.0015 

L’sr 

-.0968 

.0043 

-.0367 

.0034 

Np 

.0848 

.0019 

.0397 

.0020 

Nr 

-.3068 

.0011 

-.1925 

.0008 


-.3870 

.0005 

-.3009 

.0006 

N5a 

.0254 

.0018 

.0130 

-.0205 

.0024 

.0022 


-.2661 

.0061 

-.2555 

.0048 


-.0371 

.0008 

-.0167 

.0012 




TABLE V.- ESTIMATES OF AERODYNAMIC PARAMETERS OBTAINED 
FOR STUDY FLIGHT CONDITIONS - Concluded 


(1) a = 30°; case (b); h = 15 240 m (50 000 ft) 


Parameter 

Unmodeled 

^ effects 

Modeled 

/3 effects 

Estimated 

value 

Deviation 

Estimated 

value 

Deviation 

Lp 

-0.9764 

0.0011 

-1.0214 

0.0011 

l; 

1.9202 

.0005 

2.1204 

.0004 


-3.4293 

.0005 

-3.1528 

.0004 




-.1942 

.0022 


-1.0326 

.0016 

-1.0431 

.0015 

L’sr 

-.0615 

.0052 

-.0401 

.0049 


-.0370 

.0028 

.0427 

.0020 

n; 

-.0747 

.0014 

-.1964 

.0009 

^/3 

.1194 

.0009 

.0787 

.0009 




.0119 

.0099 

“a 

-.0556 

.0021 

-.0173 

.0018 

n’s 

L»r 

-.2360 

.0089 

-.2555 

.0063 

Y* 

^i3 

-.1487 

1 

.0009 

-.0156 

.0008 
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TABLE VI.- COMPARISON OF ACTUAL STABILITY CHARACTERISTICS 
WITH THOSE COMPUTED FROM ESTIMATION RESULTS 
FOR NON-/3-MODELING CASE 


Test 

run 

Actual characteristics 

Non-^-modeling characteristics 

Aperiodic mode 

Oscillatory mode 

Aperiodic mode 

Oscillatory mode 

t^/ 2 , sec 

^1/2^ S6C 

P, sec 

SGC 

S6C 

S0C 

P, sec 


1 

14.9 

0.429 

4.54 

1.41 

23.3 

0.365 

4.69 

1.42 

2 

5.39 

.915 

5.40 

.891 

4.47 

.588 

5.62 

.90 

3(a) 

-3.82 

.254 

53.9 

2.37 

2.66 

.342 

15.7 

1.75 

3(b) 

6.75 

1.50 

9.63 

.540 

13.9 

.855 

4.45 

.619 


39 



X 



40 


SIMULATION OF ESTIMATION OF AERODYNAMIC 

FLIGHT DATA PARAMETERS 
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Simulated flight data 
Computed using estimation results 
















O Actual value 

□ Estimated with modeled p effects 

■ Estimated without p effects 

O Combined derivatives c 
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Flagged symbols denote case (b) at a = 3(P 




Figure 9.- Comparison of estimated aerodynamic derivatives 
with actual values, h = 0. 
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Flagged symbols denote case (b) at a = 3(P 
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Figure 9.- Continued. 
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Figure 9.- Concluded. 
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Figure 10.- Continued. 
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□ Estimated with modeled p effects 

■ Estimated without p effects 

Unfiagged symbols denote case (a) at a = 3(P 
Flagged symbols denote case (b) at a = 3(P 




Figure 11.- Comparison of estimated aerodynamic derivatives 
with actual values, h = 15 240 m (50 000 ft). 
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Figure 11.- Continued. 
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Figure 11.- Concluded. 


Time , sec 


Figure 12.- Comparison of actual aileron doublet response with response computed by 
using the results from the unmodeled ^ effects estimation technique, a = 30*^; 
case (a); sea level. 
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